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The fidelity decay in a microwave billiard is considered, where the coupling to an attached antenna is varied. 
The resulting quantity, coupling fidelity, is experimentally studied for three different terminators of the varied 
antenna: a hard wall reflection, an open wall reflection, and a 50 Q. load, corresponding to a totally open chan- 
nel. The model description in terms of an effective Hamiltonian with a complex coupling constant is given. 
Quantitative agreement is found with the theory obtained from a modified VWZ approach [Verbaarschot et al, 
Phys. Rep. 129, 367 (1985)]. 

PACS numbers: 03.65.Nk, 05.45.Mt, 42.25.Bs, 03.65.Yz 



I. INTRODUCTION 

Fidelity is a standard benchmark in quantum information, 
and plays a relevant role in discussions on quantum chaos [ 1 ]. 
The corresponding fidelity amplitude can be interpreted as the 
overlap between two wave functions obtained from the prop- 
agation of the same initial state with two different time evolu- 
tions or, alternatively, as the overlap of the initial state with it- 
self after being propagated forward in time with one evolution 
and backward in time with the other. In the latter case one of- 
ten speaks of Loschmidt echo. Fidelity contains information 
on both eigenfunctions and spectra of the original and per- 
turbed systems in a non-trivial way. One can show, however, 
that there exists a profound relation [2] between fidelity decay 
and purely spectral universal parametric correlations OHa] in 
chaotic and disordered systems. The connection holds in quite 
general settings [7]. Efforts to measure fidelity are therefore 
very important. 

There was an early proposal (without the name fidelity) in 
quantum optics [8]. Along these lines the perturbation of a 
kicked rotor was discussed in great detail |9[. A realization 
of that idea is not available today, although an experiment of 
this type, but with a more complicated process, was conducted 
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Experiments with microwave cavities or elastic bodies 
seem to provide good options to study the decay of fidelity 
ITTll . but a difficulty arises. Fidelity implies an integration 
over the entire space. In two-dimensional microwave billiards 
the antenna always represents a perturbation, and thus moving 
the antenna defeats the purpose of a fidelity measurement, as 
the wave-function taken at any point is that of a slightly dif- 
ferent system. In contrast to wave function measurements, in 
fidelity experiments we are precisely interested in such differ- 
ences, and thus wave functions measured with moveable an- 
tennas lfl2l - tl4ll or a moveable perturbation body fT^ - Tnll ar e 
not appropriate. In elastic experiments on solid blocks 11181 - 
l20ll or three-dimensional (3D) microwave billiards the wave 
function inside the volume seems to be inaccessible anyway 
i2lll22ll . This leads to the development of the concept of scat- 
tering fidelity Hill which tests the sensitivity of S-matrix ele- 



ments to perturbations. This is also of intrinsic interest since 
the scattering matrix may be considered as the basic building 
block at least in the case of quantum theory 11231, 12411 ■ 

In former studies the scattering fidelity has been investi- 
gated in chaotic microwave billiards by considering a pertur- 
bation of the billiard interior. It can be shown that in such a 
case the random character of wave functions causes the scat- 
tering fidelity to represent the usual fidelity, provided that ap- 
propriate averaging is taken Ull l25ll . Scars and parabolic 
manifolds will obviously change that correspondence, but 
their effect can be avoided in experiment. Specifically, two 
different types of interior perturbations were experimentally 
studied. In the first set of experiments a billiard wall was 
shifted, realizing the so-called global perturbation Bill 12611 . 
meaning that there is a total rearrangement of both spec- 
trum and eigenfunctions already for moderate perturbation 
strengths. Good agreement with prediction from random ma- 
trix theory (RMT), expecting Gaussian or exponential decay 
depending on perturbation strength, was found. In the second 
experiment a small scatterer was shifted inside the billiard, the 
wave function being influenced only locally 12511 . Using the 
random plane wave conjecture, an algebraic decay was pre- 
dicted and confirmed experimentally. 

Actually, any measurement opens the system. Coupling 
to the continuum changes drastically the system properties 
by converting discrete energy levels into unstable resonance 
states. The latter reveal rich dynamics when the coupling 
strength to the scattering channels is varied 112711 . see also 02811 
for relevant microwave studies. Since the time evolution op- 
erator is subunitary in this case, there appears the leakage of 
the norm inside the scattering system 112911 . This decay is fully 
controlled by the degree of system openness and may also be 
considered as a remote analog of fidelity decay for open sys- 
tems. In the framework of the scattering fidelity coupling to 
the continuum can be taken into account naturally. 

It seems therefore attractive to study the sensitivity of S'- 
matrix elements to perturbations in the coupling between the 
scattering system and decay channels. This will be the central 
purpose of the present paper. Experimentally, we realize the 
system by a flat microwave billiard with two attached antennas 
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FIG. 1. Geometry of the chaotic Sinai billiard, length / = 472 mm, 
width w — 200 mm and a quarter-circle of radius r — 70 mm where 
an antenna with different terminations may be introduced at position 
c. a denotes the measuring antenna. The additional elements were 
inserted to reduce the influence of bouncing balls. 



Kit) being part of our microwave equipment. For case (a) the 
reflection amplitude S cc was also measured. From this mea- 
surement the coupling strength of antenna c can be obtained, 
see Eq. © below. For all four cases we measured 18 differ- 
ent realizations by rotating an ellipse (see Fig. [U to perform 
ensemble averages. 

An alternative to the coupling of an antenna with variable 
end is an open wave guide whose coupling to the billiard can 
be varied by a variable slit. It showed up that, contrary to in- 
tuition, for this setup the main effect of the variation of the 
slit does not correspond to a change of the coupling to the 
outside, but to a distortion of the wave functions in the bil- 
liard, thus corresponding more to the case of a local scattering 
fidelity 12511 . This system is discussed in AppendixlAl 



and measure the reflection in one antenna while modifying the 
coupling in another, see Sec.HJfor details on the experimental 
setup. Section|Ill]presents a theoretical consideration based on 
RMT and the effective Hamiltonian approach. In Sec.lIVIwe 
discuss in detail the experimental results and compare them 
with the theory. Our main findings are then summarized in 
the concluding Sec.lVl 

II. EXPERIMENT 

The basic principles of billiard experiments with mi- 
crowave cavities as a paradigm of quantum chaos research 
are described in detail in [30]. Therefore, we concentrate on 
the aspects of relevance to the present study. Reflection and 
transmission measurements have been performed in a flat res- 
onator, with top and bottom plate parallel to each other. The 
cavity can be considered as two-dimensional for frequencies 
v < Cmax = c/ (2ft), where h = 8 mm is the height of the 
resonator. 

The setup, as illustrated on Fig. Q] is based on a quar- 
ter Sinai shaped billiard. Additional elements were inserted 
into the billiard to reduce the influence of bouncing-ball res- 
onances. The classical dynamics for the chosen geometry of 
the billiard is dominantly chaotic. At position a one antenna 
is fixed and connected to an Agilent 8720ES vector network 
analyzer (VNA), which was used for measurements in a fre- 
quency range from 2 to 18 GHz with a resolution of 0.1 MHz. 
We measured the reflection S-matrix element S aa first for 
the unperturbed system, which corresponds to the situation, 
where no additional antenna is inserted at position c. Then we 
perturbed the system by inserting another antenna at position 
c which was terminated consecutively in three different ways: 

(a) connection to the VNA (total absorption), 

(b) standard open (open end reflection), 

(c) standard short (hard wall reflection), 

and again measured the corresponding reflection at antenna a 
for each case. The connection of antenna c to the VNA cor- 
responds to a termination of antenna c with a 50 load. The 
terminators for the cases (b) and (c) have been taken from the 
standard calibration kit (Agilent 85052C Precision Calibration 



III. THEORY 

A. Generalized VWZ approach to fidelity 

The general case of M scattering channels connected to N 
levels of the closed cavity can be described in terms of the 
following effective non-Hermitian Hamiltonian 

M 

H cS = H-iY,^kV k Vl. (1) 

fc=i 

Here, the internal Hamiltonian H of the closed system is rep- 
resented by a Hermitian N x N matrix, whereas Vk are M 
vectors of length N containing the information on the cou- 
pling of the levels to the continuum. The Vk are assumed to 
be normalized to one, V^Vk = 1, and Xk is the coupling con- 
stant of channel k. 

Such an approach was initially developed in nuclear physics 
i3ll - [33ll and since then has been successfully applied to study 
various aspects of open systems, including wave billiards 
II30L 134 - 13611 . Usually, the phenomenological coupling con- 
stants Afc are considered as real numbers which enter the final 
expressions via the so-called transmission coefficients. How- 
ever, in the present case of the antenna variation one has to 
consider the coupling to the variable antenna, A c , as a com- 
plex number, see discussion in Sec. lIIICl below. For the sake 
of generality, we will treat all A& as complex numbers with 
the only constraint on their real parts Re(Afc) > 0, due to 
the causality condition on the S'-matrix. We note that quite a 
similar problem of nonzero Im(Afc) arises in shell-model cal- 
culations due to the principle value term of the self-energy 
operator, cf. 13111 and 13211 . This requires proper modification 
of the theory which we briefly outline below. 

According to the general scattering formalism j3lll32Tl . the 
resonance part of the S'-matrix at the scattering energy E can 
be expressed in terms of H e g as follows: 

Sab (E) = S ab - 2i N /Re(A a )Re(A h ) V? = * V b . (2) 

Being interested in a reflection amplitude in channel a, it is 
possible, following ll36ll . to obtain another representation for 
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an arbitrary diagonal element S aa - To this end, it is conve- 
nient first to single out the contribution to H c r due to channel 
a by writing H c ff = H® s — iX a V a VJ, and then treat V a as a 
rank 1 perturbation to the term = H — i J2 k ^ a XkVkV k .. 
The upper index "a" for H£ s denotes that the contributions of 
all channels save the given one, a, are included. In the next 
step we expand (E — H^)^ 1 into a power series with re- 
spect to (E - i^ c a ff )" 1 V a V r a t and, after summing up the result- 
ing geometric series, obtain the following general relationship 
(Dyson's equation) for the corresponding resolvents IT3311 : 



E-H, 



cir 



E 



TTCl 
n Off 



E - H a . t 



ir 



l + iKVl 



1 



-0 
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H a cS 



(3) 
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This identity, being substituted in Eq. ©, yields the following 
expression for the reflection amplitude S aa : 

1 



Saa(E) = 



i-^K 1 — 



err 



1 + iXaVl 



1 



(4) 



e - m 



-v a 



ir 



Representation (0]l is very convenient to perform statisti- 
cal averaging. Adopting RMT approach to model intrinsic 
chaos, we take H from the Gaussian orthogonal ensemble 
(GOE) of N x N random real symmetric matrices which is 
the appropriate choice for the systems with preserved time- 
reversal symmetry Oil l37tl . The quantities V k are considered 
as fixed real iV-dimensional vectors of unit length. They are 
also supposed to be mutually orthogonal that ensures the ab- 
sence of the direct "fast" processes (which could be due to a 
nondiagonal part of the average S-matrix) 03211 . In the limit 
N — > oo, for finite M, the leading term for the average value 



of the resolvent is well known to be 

l(E/2) 



iy/l — E 2 /4\8 nm , where the imaginary part ac- 
counts for the famous Wigner's semicircle law. This implies 
the following result (valid up to the terms of the order of 
M/N) for the average S'-matrix JHHl, 

(Sab) = l^Sab (5) 



1 + Aa 

Here, we set E = as usual. As a result, the transmission 
coefficient takes the following form: 



T a = l-\(S aa )\ 2 = 



4Re(A a ) 
ll + AJ 2 



1,...,M, (6) 



which is in agreement with the result of Ref. M32H obtained 
by a supersymmetry calculation. It is worth noting that in the 
case of real A one gets T — n^Xp e IP> !]■ ^ n tne case 
of purely imaginary A corresponding to perfect reflection, the 
channel is closed, T = 0. 



B. Coupling fidelity 

We now proceed with the discussion of the scattering fi- 
delity. Its amplitude is defined in terms of the S'-matrix el- 



ements Fourier transformed into the time domain as follows 



fab(t) 



(s ab (t)s'; b (t)) 



(Sab(t)s ab (t))(su(t)s>: b (t)) 



(7) 



The prime indicates a change of the effective Hamiltonian of 
the original system after a small perturbation for the backward 
time evolution. Definition ^ guarantees that f a b(0) = 1. 
Furthermore, an overall decay of the correlation functions due 
to absorption drops out, provided the decay is the same for the 
parametric cross-correlation functions in the nominator and 
the autocorrelation functions in the denominator full . The 
scattering fidelity itself is 



F ab (t) = \f ab (t)\' 



(8) 



Note that in contrast to the original definition of the scattering 
fidelity, we allow for a change in the channel vectors as well. 

As it is explained above, the original idea is to change only 
the complex coupling strength A c to one channel c, while the 
measuring is done on one or two different channels a,b ^ c. 
We denote the resulting scattering fidelity by coupling fidelity . 
We present below an exact RMT prediction for this quantity. 

The starting point is to apply the convolution theorem for 
Fourier transforms to Eq. (|7]) and relate it to the parametric 
cross-correlation function C\S ao , S'* b ](t) of the S-matrix ele- 
ments in the time domain 138], 



(s ab (t)s , ; b (t))=c[s ab ,s a i}(t). 



(9) 



We denote the coupling constant for the variable antenna c 
in forward and backward time evolution by A and A', respec- 
tively. (We omit the lower index "c" henceforth.) In the case 
of unchanged coupling, A = A', the autocorrelation func- 
tion C[S ab , S* b ] (t) is real and its exact expression is obtained 
from Verbaarschot-Weidenmiiller-Zirnbauer (VWZ) integral 
ll32ll and is given by 

C[S ab ,S* ab }(t) = 6 ab T 2 {l-T a )J a (t) + (l + 6 ab )T a T b P ab {t). 

(10) 



It is convenient to use the parametrization of Ref. 0911 to 
write down the explicit expressions for the functions J a (t) and 

Pab(t), as 



J ait) = 41 



r + T a x 



t - r 



and 



Pab{t) - 21 



l + 2T a r + T 2 x l-T a {t-r) 



T a T b x 2 + d ab {r)x + (2r + l)r 



(11) 



(1 + 2T a r + T 2 x)(l + 2T b r + T^x) 
(t-r){r + l- t) 
\l-T a (t-r)][l-T b {t-r)] 



(12) 



where 



2r + 1 
2u + 1 



u 2 , d ab (r) = T a T b + (T a + T b ){r + 1) - 1 
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and the shorthand X stands for the integral, 



TV--] = 



dr I du 

max(0,t-l) JO (2 U 
M 



(t - r)(r + 1 - t) 



W 2 



x 



l-T fe (t-r) 



2 ^l + 2T k r 



(13) 



Here and below, t denotes the dimensionless time measured 
in units of the Heisenberg time tn — 2nh/A, where A is the 
mean level spacing. 

The calculation of the correlator [Eq. ®1 in the case of 
A 7^ A' proceeds along the same lines as in 03211 . see Appendix 
iBl The result turns out to be formally given by the same VWZ 
expression (TlOb . where the transmission coefficient T c in the 
varied channel c has to be substituted by 



2(A + A'*) 



(14) 



(1 + A) (1 + A'*) ' 

while performing the integration [Eq. (H3H. The quantity Tj? s 
may be considered as an effective transmission coefficient due 
to a parametric variation of the coupling strength in the chan- 
nel c. Only if A — A', Tjr 1 becomes equal to the conventional 
transmission coefficient [Eq. In contrast to Eq. ©, T£ s 
is generally complex and also T c cff ^ 1— (S cc ) (S'* c ) . We note, 
however, that Eq. (fl4] i can be cast in the following form 



l c — 1 — o cff o c 



If ' 



where 



S'cff 



1 - A'* 
1 + A 



J cff 



1 - A 

1 + A'* 



(15) 



(16) 



These quantities might be interpreted as the (average) para- 
metric S'-matrix amplitudes in the varied channel for the for- 
ward and backward time evolution, respectively. 

The subsequent evaluation of coupling fidelity cannot be 
done analytically and will be performed numerically. 



C. Effective Hamiltonian description 

The experimental situation shall now be mapped onto the 
theory derived in the preceding subsection. Though the cal- 
culation is straightforward, and similar approaches can be 
found elsewhere [40], it is repeated here for the reader's con- 
venience. Let us start with the expression of the scattering 
matrix in terms of Wigner's reaction matrix: 



S = 



1 - iW^GW 



1 



iW^GW 



(17) 



G = (E — H) -1 is the Green's function of the closed system 
and matrix W = (W a ,W c ) contains the information on the 
coupling. As before, index "c" refers here to the antenna with 
variable coupling, and "a" to the measuring antenna. Per def- 
inition, the 5-matrix relates the amplitudes of the incoming 
(it) and outgoing (v) waves, 



S 



(18) 



A termination of antenna c is described by 



u c = rv c , r = e- {a - lv) . 



(19) 



where r contains the information on the reflection properties 
of the antenna. For reflection at an antenna with open or 
closed end we have a — (as long as the absorption in the 
antenna can be neglected). The termination of the antenna by 
a 50 fl load corresponds to a — > oo. 

Making use of Eq. ( fTTb . one can rewrite Eq. ( fT~8T > as 



iW^GW 



( U c - 




1 = 1 


f U c -V c \ 


\ u a - 


\-v a J 




[ u a -v a J 



(20) 



Substituting relation ( fT9l ) in Eq. d20i >, u c and v c can be elimi- 
nated, resulting in an equation for u a and v a , 



iWlGW a ( 



(21) 



Here, we have introduced the modified Green's function, G, 
with the following matrix element 



WlGW a = G aa - G ac 



1 



1 + i\ T G c 



-iXtGc 



(22) 



where G nm — W^GW m and At is the coupling constant of 
the "terminator," 



A, 



1 



tanh ■ 



(23) 



Equation (fJTJ has the same form as Eq. (l20l . but for the 
measuring antenna only and with the modified Green's func- 
tion. Substituting explicit expressions for matrix elements 
Gnm, we obtain in a number of elementary steps 



G = G- 



1 



1 



1 + iXrWcWjG 



(24) 



where iJ" ff = H — iXrW c W^ . Introducing the normalized 
coupling vector V = ^= W c , where \w = W\W C is a 
channel coupling strength, Hg s may be finally written as 

H^g = H-i\VV\ \ = \ T \ W - (25) 

The total coupling constant A is generally complex and takes 
into account the effects of both the channel coupling (\w) and 
the terminator (At)- The 2x2 scattering matrix [Eq. ( TT7i il for 
the measuring antenna and the antenna with variable termina- 
tor has thus been reduced to a 1 x 1 scattering matrix for the 
measuring antenna only, 



1 - iWl 



1 



a E- 



Tja 
cti 



-W n 



l+lW a 



1 



(26) 



e - m, 



■w a 



In the case of a single measurement antenna and one antenna 
with variable coupling, Eq. d26b is equivalent to Eq. (|4j. Equa- 
tions d24"li-(|2"6ii constitute the main result of this section. They 
show that the influence of the variable antenna can be taken 
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into account by an appropriate modification of the Hamilto- 
nian. 

Two special cases are of particular importance. For the ter- 
mination of the antenna with a 50 il load the outgoing wave is 
completely absorbed, corresponding to the limit a — > oo. It 
follows At = tanh oo = 1 and 

ff c a ff = H - iX w VV^ . (27) 

In this case the coupling is purely imaginary. For the two 
cases, where the antenna is terminated by a hard wall or an 
open reflecting end, we may assume a = 0, resulting in At = 

tanli(i</j/2) = ztan^/2, and 

H« s = H + tan (|) X W VV ] . (28) 

In this case the coupling is purely real, and the antenna does 
not correspond any longer to an open channel but to a scatter- 
ing center only. This is true, as long as the absorption in the 
antenna can really be neglected. This becomes questionable, 
as soon as <p approaches ir, corresponding to the excitation of 
a resonance within the antenna. For this singular situation the 
perturbative treatment of the antenna coupling applied in the 
derivation looses its justification. 

The value of ip depends on the length of the antenna in units 
of the wave length and thus on frequency. But independently 
of frequency the difference of the phase shift ip for the reflec- 
tion at the open end (oe) and the hard wall (hw), respectively, 
is always ir. A phase difference of ir means a replacement 
of the tangent by the cotangent in Eq. ( 1281 . i. e. the coupling 
constants At for the two situations are related via 

A T ,hw = 1/Atoc (29) 

With the above introduced total coupling constant A = 
AtAvk> this may be alternatively be written as 

AhwAoo = = ^50fi (30) 

since Xw is the coupling constant for the 50 £1 load, see 
Eq. ( 1771 ). Ahw and A oc denote the total coupling constants for 
the hard-wall and the open-end reflections. These relations 
allow for explicit tests of the theory. 

IV. RESULTS AND DISCUSSION 

In this section we want to discuss the experimental and the- 
oretical results for the coupling fidelity decay under the per- 
turbations (a)-(c) described in Sec. [II] For all results below the 
system without the varied antenna, corresponding to A = 0, 
is chosen as the reference, whereas for the perturbed system 
the coupling constant is A' = Asorj, A O0 , or Ahw, depending 
on the terminator. 

We start with a plot of the complex valued fidelity ampli- 
tude for one frequency range, see Fig. [2] The solid lines show 
the experimental results, derived from the Fourier transform 
of the measured S aa (v) and S' aa (v) via relation (0 for the sit- 
uations (a) 50 £1 load (black), (b) open-end reflection (green, 
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FIG. 2. (color online) Real part fnit) and imaginary part fi(t) of 
the fidelity amplitude for three types of perturbation: Ason (black); 
Ahw (orange, light gray); A oc (green, dark gray) in frequency range 
8.0 — 8.5 GHz. The time is given in units of the Heisenberg time 
tn = 2-irh/A, where A is the mean level spacing. Solid lines show 
the experimental results. The theoretical curves are dotted for exper- 
imental parameter (available only for the 50 load case) and dashed 
for fitting parameter. The corresponding parameter and the transmis- 
sion coefficient for the measuring antenna a are listed in Tab. U The 
black dotted curve is nearly indistinguishable from the dashed one. 

dark gray) and (c) hard-wall reflection (orange, light gray). 
For the case (a) we are able to calculate the corresponding 
theoretical curve (black dotted line) without any free parame- 
ter, since the coupling constant Xw can be determined directly 
from the additional reflection measurement at antenna c (see 
Sec. HJi using relation ©■ According to Eq. d27b we expect 
a purely imaginary coupling with X50Q. For our antenna Xw 
varies from 0. 1 to 0.4 in a range from 6 to 10 GHz. Using the 
experimentally determined Xw, one gets already very good 
agreement between experiment and theory for the 50 £1 load 
without any fit. A fit of A^ to the experimental curves only 
marginally improves the correspondence. 

For the perturbations with reflecting ends the systems are 
closed. The correct description for this situations is Eq. d28l . 
With the total coupling constants A oc and Ahw as a free fitting 
parameter we find again an agreement between experiment 
and theory (dashed lines), for both, real and imaginary part. 
As one would expect from theory for the case of reflecting 
ends we see a significant imaginary part of the fidelity ampli- 
tude fi(t) (green and orange lines), whereas in the case of an 
absorbing end (black line) the imaginary part is nearly zero. 

It is convenient to continue our discussion in terms of the 
fidelity (not its amplitude) introduced in Eq. (0. In Fig. [3] we 
present the experimental and theoretical fidelity results for the 
situations 50 il load (black lines), open-end (green lines), and 
hard-wall (orange lines) as perturbation, for three frequency 
ranges. In case of the closed channels the total phase shift 
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FIG. 3. (color online) Fidelity decay |/(i)| 2 for three types of per- 
turbation: Xson (black); Ah w (orange, light gray); A oc (green, dark 
gray) in three different frequency ranges: (i) 7.2 — 7.7 GHz; (ii) 
8.0 - 8.5 GHz; (iii) 8.7 - 9.2 GHz. Solid lines show the experi- 
mental results and the theoretical curves are dotted for experimental 
parameter (available only for the 50 Q load case) and dashed for fit- 
ting parameter. The corresponding parameter and the transmission 
coefficient for the measuring antenna a are listed in Tab. [J 



<f(y) increases monotonically with frequency. Thus A oc and 
Ahw are oscillating in counter phase, see Eq. (f28t . This in- 
duces a corresponding oscillation in the strength of the fidelity 
decay as seen in Fig. [3] For a more quantitative discussion we 
compare the experimentally determined fidelity decay (solid 
lines) to the theoretical curves (dotted and dashed lines). First 
of all, one sees that fitting (dashed lines) the experimental re- 
sults again works well for all cases. Focusing on the 5017 
load case (black lines) also the theoretical results without free 
parameter (dotted lines) show good agreement with the ex- 
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\ cxp 
A 50f2 




A oc 


Ahw 


\w 


T a 


0) 


7.2 - 7.7 


0.19 


0.37 


0.65 % 


— 0.04 % 


0.16 


0.19 


(ii) 


8.0-8.5 


0.21 


0.20 


0.19z 


— 0.23 % 


0.21 


0.22 


(iii) 


8.7-9.2 


0.24 


0.21 


0.05 1 


— 0.83 1 


0.20 


0.34 



TABLE I. Coupling constants in the three different frequency ranges 
(i)-(iii). According to Eq. BOb . AgoQ and Agon should be compared 
to Xw = V-^ocAhw, see main text for discussion. 

periment for the frequency ranges plotted in Figs. [3] (ii) and 
(iii). Only the plot in Fig. [3] (i) shows significant deviation 
between the theoretical result without free parameter and the 
experimental curve. We want to stress that the case shown 
here is the worst among all the investigated frequency ranges. 
Here the experimental fidelity amplitude /son shows a signif- 
icant imaginary part. Thus the imaginary part of A500 is not 
zero, i.e. the 50 £1 terminator does not correspond to perfect 
absorption, and Eq. ( f28T > does not hold. This might be due to 
an antenna resonance, leading to an increased reflection from 
the channel c. 

Finally we perform a check on the coupling constants based 
on Eq. d30b . Accordingly, the square root of the product 
of the coupling constants for open-end and hard-wall reflec- 
tion should give Agon. Table U shows that for the frequency 
ranges (ii) and (iii) there is indeed good agreement between 
VAocAhw, AfJ n and A^. In the case (i) VAocA hw agrees 
quite good with the experimental parameter A^^, but the fit- 
ting parameter is much larger. This deviation reconfirms our 
arguments presented in the above discussion of the fidelity 
plot shown on Fig. |3ji). 



V. CONCLUSIONS 

In this work, we have studied the influence of the coupling 
to the continuum on the decay of fidelity. This complements 
previous experiments of our group, where the fidelity decay 
under the influence of various types of geometrical perturba- 
tions was studied ITlll25ll26ll4lll but for closed systems exclu- 
sively. To get rid of an overall absorption we used the concept 
of scattering fidelity introduced by us previously Hill , defined 
as the parametric cross-correlation function of S-matrix ele- 
ments normalized to the corresponding autocorrelation func- 
tion. 

On the theoretical side we have developed a model descrip- 
tion of the fidelity decay in terms of a modified VWZ ap- 
proach. The parametric cross-correlation function of 5-matrix 
elements for two different A ^ A' can be reduced to an auto- 
correlation function with a complex effective transmission co- 
efficient [Eq. dl4H . thus expressing coupling fidelity in terms 
of a modified VWZ integral. This theory holds for an arbitrary 
number of channels and describes the experimental coupling 
fidelity results well. It would be interesting to investigate, 
whether it is possible to relate coupling fidelity obtained via 
the VWZ ansatz to an approximation using the random plane 
wave conjecture, as used to explain the local fidelity B25I1 . 

We have found two additional important result. First, a 
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smooth variation of the coupling, e. g. by varying the coupling 
to an attached wave guide will not easily yield the information 
about the effect of coupling to the continuum on the scatter- 
ing fidelity. Each geometric variation will give rise to both 
a change of coupling and internal scattering properties, thus 
screening the purely external effect, as is also discussed in 
AppendixlAl 

Second, we have included closed channels within the de- 
scription of VWZ. The speed of the fidelity decay for the 
open-end and hard-wall reflection oscillates with frequency 
due to the corresponding variation of the phase with fre- 
quency. An important relation (f30b between the coupling con- 
stants of the antenna terminations has been established, en- 
abling us to connect the results found for the closed channel 
(A oc , Ahw) to those for the open channel (A500). In all cases 
the fidelity decay for at least one of the reflecting antennas is 
faster than for the open channel, showing the strong influence 
of the imaginary part on the coupling constant A. 
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Appendix A: Experiment with attached wave guide with 
variable coupling 

In our first approach we used the setup shown in Fig. [4] 
where the opening of the variable slit plays the role of the fi- 
delity parameter. The setup is based on a quarter Sinai shaped 
billiard with length I — 342 mm, width w = 237 mm and 
a quarter-circle of radius r = 70 mm, and an attached chan- 
nel. The channel has a total length l c — 243 mm and a width 
w c = 16 mm. At position a and c two antennas were fixed and 
connected to the VNA. The complete S-matrix was measured 
in a frequency range from 9.5 to 18.0 GHz with a resolution of 
0.1 MHz, where the wave guide only supports a single prop- 
agating mode, i. e. it acts as a single channel. The perturba- 
tion of the system was achieved by opening the channel from 
d = — 16 mm in steps of 0.1mm using a slit diaphragm 
at the point of attachment. An ellipse insert with semiaxis 
a = 70 mm and b = 40 mm was rotated to get an ensemble of 
20 different systems for averaging. Additional elements were 
inserted into the billiard to avoid bouncing-ball resonances. 
The wave guide was terminated by a perfect absorber, which 
according to Eq. (l27l > should correspond to a purely imagi- 




100 mm 



FIG. 4. Geometry of the billiard with attached wave guide. 

nary coupling. As before the coupling constant Xw could be 
determined directly from a reflection measurement at antenna 
c. Aw could be varied from Xw — (no coupling) to Xw = 1 
(perfect coupling) by increasing the opening d of the slit. 

In Fig. |5] the coupling fidelity decay is shown for two dif- 
ferent perturbation strengths. The solid lines show the exper- 
imental results. With the formulas derived in Sec. [TIT] we cal- 
culated the expected theoretical fidelity decay assuming that 
the channel is totally open, i. e. the coupling is purely imag- 
inary (dotted lines). There is obviously no agreement. This 
shows that something is wrong in the argumentation. For a 
further check we removed the absorbing end and the antenna 
in the channel and replaced it by a reflecting end thus closing 
the system. We did not find any noticeable difference to the 
case with the absorbing end and the antenna in the channel ex- 
perimentally. So there is only one explanation: by far the ma- 
jor part of the wave is reflected directly at the slit, and only a 
minor part really penetrates into the channel. This means that 
the coupling is not imaginary but mainly real (up to perhaps 
a minor imaginary contribution), and we should use Eq. ( 1281 




0.0 1 1 , 1 , 1 , , 1 , 1 , 1 , 1 1 , 1 , , 

0.0 0.5 1.0 1.5 2.0 

t 

FIG. 5. (color online) Experimental coupling fidelity \f{t)\ 2 (solid 
lines) and theoretical results for the experimental parameter Xw (dot- 
ted line) and the fit parameter Aflt (dashed line) for two openings 
d — 6.5 mm with Xw ~ 0.05 and Aflt = — 0.18z (black), and 
d = 11.2 mm with Aw = 0.52 and Aflt = — 0.55i (orange, light 
gray). The frequency window of the Fourier transform of the mea- 
sured Sa,a,(v) was 13 to 14 GHz and the transmission coefficient for 
antenna a was T a = 0.95. 
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instead of Eq. (|27| | for the interpretation of our results. The 
dashed lines in Fig. [5] show the resulting theoretical curves 
with Afit as a free parameter according to the definition pre- 
ceding Eq. d28l >. Now a perfect agreement between experi- 
ment and theory is found. 

As a resume we can state that the variable slit works essen- 
tially as a scattering center leading to partial masking of the 
change of coupling by the change of scattering properties in 
the fidelity decay. 



Appendix B: Derivation of EQ. fl4t and qualitative discussion 

The calculation of Eq. (0 proceeds along the same lines as 
in l32ll . hence we indicate below only the main steps and es- 
sential differences. First, we make use of the representation 
of resolvents and thus S'-matrix elements [Eq. d2J] in terms of 
Gaussian integrals over auxiliary "supervectors" consisting of 
both commuting and anticommuting (Grassmann) variables. 
This allows us to perform statistical averaging over GOE ex- 
actly. Then in the RMT limit N — >• 00, the remaining integra- 
tion over the auxiliary field can be done in the saddle-point ap- 
proximation. It turns out that there exists a nontrivial saddle- 
point manifold ll42ll over which one has to integrate exactly. 
As a result, the two-point correlation function of the 5-matrix 
elements in the energy domain acquires the form of a certain 
expectation value in field theory (nonlinear supersymmetric 
cr-model), ((• ■•)) = J A[a G ]e c ^ F M {- ■ ■)> cf. Eq. (7.13) of 
Ref. 13211 . In the notations of this paper, the effective La- 
grangian reads C(e) — jNetig(aGL), with e being the en- 
ergy difference. Definitions of the supertrace, trg, as well as 
of the supermatrices <jq and L can also be found there (see 
ll43ll for a general reference). The pre-exponential terms omit- 
ted above depend on the coupling constants in the channels a 
and b (a,b ^ c), being thus the same as considered in l32ll . 
They finally correspond to the expressions appearing explic- 
itly in Eqs. ( fTTT i and ( fT2l . At last, the so-called channel factor 
J~m accounts for the coupling to all the channels. It is the 
term that requires modification due to both generally complex 
and varied coupling constants. In the [1,2] block notation (the 
"advanced-retarded" ordering of supermatrix elements), Tm 
reads 



M 



M ( l 

= n ° x p -2 trgin 
fc=i 1 




a G L 



where = A', for all channels save the varied one, k 7^ c. 
By employing the "angular" parametrization of ac in terms 
of the matrices t\2 and £21, the subsequent evaluation of J 7 ^/ 
goes in parallel with Sec. 7 of Ref. 13211 . with the final result 
being 



M 



M 

n ex p 

fe=i 



-7; trg ln(l 



Tfu,U 



(Bl) 



This is just a usual formula for the channel factor in the VWZ 
theory except for the effective transmission coefficient T° in 
the channel c that is now given by expression dT4b (we note 
that Tf = T k if k ^ c) & Performing finally the Fourier 
transform, the two-point correlation function in the time do- 
main takes the form of Eqs. (ITOt— (fl~3l>. with the above mod- 
ification in the channel factor corresponding explicitly to the 
second line of Eq. ( [PTl i. 

Although the subsequent evaluation cannot be made analyt- 
ically and has to be done numerically, it is still useful to make 
some qualitative analysis. To this end we note that P a b(t) 
and J a (t) are quite similar in structure to the "norm leak- 
age" decay function [29] and the form factor of the Wigner 
time delays B45I1 . Following the analysis performed there 
(see also 03511 ). one notices that the time dependence in ques- 
tion is mainly due to the channel factor [Eq. dBlH . In the 
time domain, its typical behavior is ~ Ofe=i(l + f^fc*) ~^ 2 , 
where (3 = 1 is for the present case of time-invariant sys- 
tems whereas (3 — 2 is for the case of broken time invari- 
ance (GUE). The case of the coupling variation in the channel 
c amounts then to replacing in this expression T c with T° ff . 
This suggests the following heuristic formula for the coupling 
fidelity 



(l + 2T c f//?)(l + 23£/0) 



\l + 2Tft/f3\ 2 



0/2 



(B2) 



For the parameters of A and A' found in the experiment only 
deviations on the % level were found while making fit to sur- 
mise dB2t . We stress, however, that there is no control on ap- 
proximations involved to derive this expression. One should 
generally expect that surmise (1B2I ) coincides closely with the 
exact result at small times (when it is given by an exponen- 
tial dependence), while the exact asymptotic behavior at large 
times is reproduced up to a factor of the order of unity (as was 
indeed confirmed numerically). Therefore, we have used the 
exact supersymmetry result for all the figures and analysis of 
the main text. 
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